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Forward-backward multiplicity correlation strengths have been measured with the STAR detector 
for Au+Au and p+p collisions at ^/SNN = 200 GeV. Strong short and long range correlations (LRC) 
are seen in central Au+Au collisions. The magnitude of these correlations decrease with decreasing 
centrality until only short range correlations axe observed in peripheral Au+Au collisions. Both the 
Dual Parton Model (DPM) and the Color Glass Condensate (CGC) predict the existence of the 
long range correlations. In the DPM the fluctuation in the number of elementary (parton) inelastic 
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collisions produces the LRC. In the CGC longitudinal color flux tubes generate the LRC. The data is 
in qualitative agreement with the predictions from the DPM and indicates the presence of multiple 
parton interactions. 

PACS numbers: 25.75.Gz 



The study of correlations among particles produced 
in different rapidity regions may provide an understand- 
ing of the elementary (partonic) interactions which lead 
to hadronization. Many experiments show strong short- 
range correlations (SRC) over a region of ~ ± 1 units in 
rapidity [1, 2]. In high-energy nucleon-nucleon collisions 
(\/s ^ 100 GeV) the nonsingly diffractive inelastic cross 
section increases significantly with energy and the mag- 
nitude of the long-range forward-backward multiplicity 
correlations (LRC) increases with the energy [2]. These 
effects can be understood in terms of multiparton inter- 
actions [1]. 

In high energy nucleus-nucleus collisions, it has been 
predicted that multiple parton interactions would pro- 
duce larger long-range forward-backward multiplicity 
correlations that extend beyond ± 1 units in rapidity, 
compared to hadron-hadron scattering at the same en- 
ergy 0, S [1]. The model based on multipomeron ex- 
changes 

(Dual Parton Model) predicts the existence of long range 
correlations [1, Q • In the Color Glass Condensate 
(CGC) picture of particle production the correlations of 
the particles created at early stages of the collisions can 
spread over large rapidity intervals, unlike the particles 
produced at later stages. Thus the measurement of the 
long range rapidity correlations of the produced particle 
multiplicities could give us some insight into the space- 
time dynamics of the early stages of the collisions [6( . 

One method to study the LRC strength is to mea- 
sure the magnitude of the forward-backward multiplic- 
ity correlation over a long range in pseudorapidity. 
Such correlations were studied in several experiments 
I>M_S H, I pi [lol [iTI and investigated theoretically 
I 111, Mill 

]lq|. In this paper we present the 
first results on the forward-backward (FB) multiplicity 
correlation strength and its range in pseudorapidity in 
heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) measured by the STAR experiment. Ear- 
lier analyses in STAR have focused on the relative cor- 
relations of charged particle pairs on the difference vari- 
ables A77 (pseudorapidity) and (azimuth). It was 
observed that the near-side peak is elongated in Arj in 
central Au-|-Au as compared to peripheral coUisions @. 
In the present work the measure of correlation strength 
as defined in Eq. (1) and the coordinate system differs 
from that of these earlier STAR measurements. The FB 
correlation strength is measured in an absolute coordi- 
nate system, where 77 = is always physically located at 
midrapidity (the collision vertex) , instead of the relative 
r] difference utilized in other 2-particle analyses. These 
differences allow the determination of the absolute mag- 
nitude of the correlation strength. 



The correlation strength is defined by the dependence 
of the average charged particle multiplicity in the back- 
ward hemisphere, (Nh), on the event multiplicity in the 
forward hemisphere, Nf, such that (Ni,) = a+bNf, where 
a is a constant and b measures the correlation strength: 

{NfN,)~{Nf){N,) _ Dl^ 

In Eq. (1), Z?^^ (covariance) and (variance) are the 
backward-forward and forward-forward dispersions, re- 
spectively [3,^. 

The data utilized for this analysis are from year 2001 
(Run II) y/sNN — 200 GeV minimum bias Au-f Au col- 
lisions (~ 2.5x10® events) at the RHIC, as measured by 
the STAR experiment [T3| • The FB correlation has been 
studied as a function of the centrality of the collision. 
The centralities studied in this analysis account for 0- 
10, 10-20, 20-30, 30-40, 40-50 and 50-80% of the total 
hadronic cross section. All primary tracks with distance 
of closest approach to the primary event vertex < 3 cm 
in the Time Projection Chamber (TPC) pseudorapidity 
range I77I < 1.0 and with transverse momentum > 0.15 
GeV/c were considered. This region was subdivided into 
bins of width rj = 0.2. The FB intervals were located 
symmetrically about midrapidity (ry = 0) with the dis- 
tance between bin centers (pseudorapidity gap A77): 0.2, 
0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, and 1.8. To avoid a bias 
in the FB correlation measurements, care was taken to 
use different pseudo-rapidity selections for the centrality 
determination which is also based on multiplicity. There- 
fore, the centrality determination for the FB correlation 
strength for A77 — 0.2, 0.4 and 0.6 is based on the multi- 
plicity in 0.5 < I77I < 1.0, while for A77 = 1.2, 1.4, 1.6 and 
1.8 the centrality is obtained from jryj < 0.5. For Ary = 
0.8 and 1.0 the sum of multiplicities from \ri\ < 0.3 and 
0.8 < I77I < 1.0 is used for the centrality determination. 
The effect of centrality region selection on FB correla- 
tion strength was also studied by narrowing the region to 
I77I < 0.3, 0.2 and 0.1 for all Arj bins. This increases the 
FB correlation strength by ^ 10-15% at the most. Since 
the pseudorapidity particle density {dN/drf) plateau at 
^smn = 200 GeV in Au-|-Au collisions extends over the 
region of interest [l^, this procedure yields consistent 
particle counts in the FB measurement intervals. An 
analysis of the data from (Run II) p+p collisions at ^/s 
— 200 GeV, was also performed on minimum bias events 
(^ 3.5x10® events) using the same track cuts as for the 
Au+Au analysis. Corrections for detector geometric ac- 
ceptance and tracking efficiency were carried out using a 
Monte Carlo event generator (HIJING) and propagating 
the simulated particles through a GEANT representation 
of the STAR detector geometry. 
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In order to eliminate (or at least reduce) the effect of 
impact parameter (centrality) fluctuations on the mea- 
surement of the FB correlation strength, each relevant 
quantity {Nf, Ni,, Nj, NfNf,) was obtained on an event- 
by-event basis as a function of the event multiplicity, Nch- 
The average uncorrected mean multiplicities {Nf)uncorr, 

{Nb)uncorr, {N'j)uncorr, and {NfNb)uncorr in each CCU- 

trality bin were calculated from a fit to the Nch depen- 
dences [19, 20]. The functional forms that were used are 
linear in Nf, Nt, and quadratic in Nj and NfNi, for all At] 
bins. Tracking efficiency and acceptance corrections were 
then applied to {Nf)uncorr, {Nb)uncorr, {N'j)uncorr, and 
{NfN},)uncorr to cach cvcut. Then the corrected values of 
(Nf), (Nb), (iV|), and (NfNb) for each event were used 
to calculate the backward-forward and forward-forward 
dispersions, D^j: and Djf, binned by centrality, and 
normalized by the total number of events in each bin. 
This method removes the dependence of the FB corre- 
lation strength on the width of the centrality bin. As 
a cross check an alternative method of centrality deter- 
mination was also carried out using the STAR Zero De- 
gree Calorimeter (ZDC) for the 0-10% centrality range 
and the results are shown in Fig. la along with the 0- 
10% most central events from the minimum bias dataset. 
Statistical errors are smaller than the data points. Sys- 
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FIG. 1: (a) FB correlation strength for 0-10% (circle) and 
ZDC based centrality ( square) (b) FB correlation strength 
for 10-20, 20-30, 30-40, 40-50 and 50-80% ( square) Au+Au 
and (c) for p+p collisions as a function of ISj] at ^Jsnn ~ 200 
GeV. The error bars represent the systematic point-to-point 
error. The boxes show the correlated systematic errors. 

tematic effects dominate the error determination. The 
systematic errors are determined by binning events ac- 
cording to the z-vertex in steps of 10 cm varying from 



-30 to 30 cm and the maximum value of the fit range (0- 
570, 0-600 and 0-630) for {N f) , (Nb) , (Nj) , and (NjNb) 
vs Nch- An additional error could arise due to finite de- 
tection efficiency in the TPC. This is estimated to be ~ 
5% for most central collisions. The overall systematic er- 
rors due to the fit range, which causes a correlated shift 
along the y-axis, are shown in figures as boxes. 

Figure [T] shows the FB correlation strength as a func- 
tion of Ar] for p+p and centrality selected Au+Au colli- 
sions along with the ZDC based centrality results. The 
results from ZDC are slightly lower as compared to the 
0-10% most central events sampled from minimum bias 
datasets using Nch- It is observed that the magnitude of 
the FB correlation strength decreases with the decrease 
in centrality. The FB correlation strength evolves from 
a nearly flat function for 0-10% to a sharply decreasing 
function with A77 for the 40-50 and 50-80% centrality 
bins, which is expected if only short range correlations 
(SRC) are present [1]. The FB correlation strength val- 
ues for 40-50 and 50-80% centrality bins at large gap 
(A77 > 1.0) have an average value near zero. The in- 
dividual b values are near zero within their systematic 
errors. Figure [1] shows that the dependence of the FB 
correlation strength with A77 is quite different in central 
Au+Au compared to p+p collisions. It is also observed 
that the FB correlation strength decreases faster in the 
peripheral (40-50 % centrality) Au+Au as compared to 
p+p collisions. This indicates that the short range cor- 
relation length is smaller in Au+Au collisions than in 
p+p. 

Figure 2 shows the dependence of the dispersions Dfj 
and D^j: on A77 for central Au+Au collisions (Fig. 2a) 
and p+p collisions (Fig. 2b). The nearly constant value 
of Dj:j with A77 represents the dispersion within the same 
77 window, which has approximately the same average 
multiplicity for all Ar] values. The behavior of D^^ is 
similar to the FB correlation strength. Thus the FB cor- 
relation variation with the size of A77 is dominated by the 
in Eq. (1). 

Short range correlations have been previously observed 
in high energy hadron-hadron collisions [l|. The shape 
of the SRC function is symmetric about midrapidity and 
has a maximum at 77 = 0. It can be parameterized as 
oc exp(— A77/A), where A is the short range correlation 
length and is found to be A ~ 1. Thus the SRC are sig- 
nificantly reduced by a separation of ~ 1.0 units of pseu- 
dorapidity [2l| . The short range correlation length is 
smaller in nucleus-nucleus collisions as compared to high 
energy hadron-hadron collisions [1, [l^ . The remaining 
portion of the correlation strength can be attributed to 
the LRC. This can be seen in Fig. lb where the magni- 
tude of the FB correlation strength is zero for A77 ~ 1 
for 40-50% centrality. In case of 0-10% Au+Au collisions 
the magnitude of FB correlation strength is 0.6, indicat- 
ing that 60% of the observed hadrons are correlated. 

The FB correlation results are compared with the pre- 
dictions of two models of A+A collisions widely used 
at RHIC energies - HIJING [22] and the Parton String 
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FIG. 2: (a) Backward-forward dispersion, -Dt/, and forward- 
forward dispersion, -D//, for 0-10% centrality as a function 
of A77 for Au+Au collisions at y^s^ = 200 GeV. (b) D^f, 
and D^f for p+p collisions at -^s — 200 GeV. The error bars 
represent the systematic point-to-point error. The boxes show 
the correlated systematic errors. 



Model (PSM) [23]. The PSM is the Monte Carlo im- 
plementation of the Dual Parton Model (PPM) |5i] or 
Quark-Gluon String Model (QGSM) concepts [Ij], con- 
sidering both soft and semihard components on a par- 
tonic level. The HIJING model is based on perturbative 
QCD processes which lead to multiple jet production and 
jet interactions in matter (2^ . Nearly 1 million min- 
imum bias Au-|-Au collisions at ■y/ijv]v=200 GeV were 
simulated for each model. The PSM events were ob- 
tained without string fusion options. We used HIJING 
version 1.383 with default options. We have also simu- 
lated 10 million p+p minimum bias events at the same 
cms energy to provide the reference for comparison with 
Au-f Au collisions. The correlation analysis was carried 
out exactly in the same manner as for the data. Both 
PSM and HIJING predictions are lower than the data 
as shown in Fig. 3 but PSM exhibits a large LRC for 
Arj > 1.0 while HIJING predicts significantly smaller 
correlations than observed in the data. In case of p+p 
collisions the HIJING prediction agrees with the data. 
The PSM does not show the decrease of b with Ar] as 
seen in the data. These trends are illustrated in Fig. 3, 
where the variation of the FB correlation strength with 
A77 is shown for Au+Au, HIJING, and PSM. The strong 
fall of b with A77 in HIJING provides some constraints 
on the contribution of impact parameter fluctuations to 
the correlation strength (Fig. 3). Recently, Hadrod-string 
dynamics (HSD) transport model [1^ and CGC [1^ have 
addressed the possible effect of impact parameter fluctu- 
ations on the correlations with different results. 

A description of the FB correlations, which qualita- 



tively agrees with the measured behavior of FB corre- 
lation strength is provided by the Dual Parton Model 
(DPM) [3]. As mentioned earlier the FB correlation 
strength is controlled by the numerator of Eq. [1] For 
the case of hadron-hadron collisions: 



{NfNt)-{Nf){Nt) = 
{k){{NofNob)-{Nof){Nob)) 
+ [i{k')-{kmNof){Nob) 



(2) 



where {Nof) and (A^ob) are the average multiplicity of 
charged particles produced in the forward and backward 
hemispheres in a single elementary inelastic collision 
The average number of elementary (parton-parton) in- 
elastic collisions is given by (fc). The first term in Eq. 2 
is the correlation between particles produced in the same 
inelastic collision, representing the SRC in rapidity. The 
second term, (k^) — (fc)^, is due to the fluctuation in the 
number of elementary inelastic collisions and is controlled 
by unitarity. This term gives rise to LRC 0, [1] . 

Recently, long range FB multiplicity correlations have 
also been discussed in the framework of the CGC /glasma 
motivated phenomenology [2l|, [23| . The glasma provides 
a QCD based description which includes many features 
of the DPM approach, in particular the longitudinal ra- 
pidity structure [11]. This model predicts the growth of 
LRC with collision centrality [2l|. It has been argued 
that the long range rapidity correlations are due to the 
fluctuations of the number of gluons and can only be 
created at early time shortly after the collision [^, [2§] . 

In summary, this is the first measurement of long-range 
FB correlation strengths in ultra relativistic nucleus- 
nucleus collisions. A large long range correlation is ob- 
served in central Au+Au collisions that vanishes for 40- 
50% centrality. Both DPM and CGC argue that the 
long range correlations are produced by multiple parton- 
parton interactions 0, i]. Multiple parton interactions 



0-1 0°/o Au+Au 200GeV 




Ar| 



FIG. 3: The FB correlation strength for 0-10% most cen- 
tral Au+Au collisions and p+p from data ( circle), HIJING 
( triangle) and PSM ( square). The error bars shown are for 
data. 
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are necessary for the formation of partonic matter. It 
remains an open question whether the DPM and CGC 
models can describe the LRC reported here and the 
near-side correlations simultaneously. Further stud- 
ies of the forward-backward correlations using identified 
baryons and mesons as well as the dependence of the 
correlations on the collision energy may be able to dis- 
tinguish between these two models. 
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